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Abstract. Observations of the precessing orbits of stars very near the massive black 
hole in the galactic center could provide measurements of the spin and quadrupole 
moment of the hole and thereby test the no-hair theorems of general relativity. Since 
the galactic center is likely to be populated by a distribution of stars and small black 
holes, their gravitational interactions will perturb the orbit of any given star. We 
estimate the effects of such perturbations using analytic orbital perturbation theory, 
and show that for a range of possible stellar distributions, and for an observed star 
sufficiently close to the black hole, the relativistic spin and quadrupole effects will be 
larger than the effects of stellar cluster perturbations. Our results arc consistent those 
from recent numerical A'^-body simulations by Merritt et at. 
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1. Introduction and summary 

The center of our galaxy has become an active arena for studying possible tests of 
general relativity (GR) in the strong-field regime, because of the near certainty that 
it harbors a 4 million solar-mass black hole, colloquially denoted SgrA* (see [U E] for 
reviews). Numerous authors have studied the observability of relativistic effects in 
the vicinity of the black hole, in the observable motion and behavior of orbiting stars 
[31 m El El [71 m El [ini E], in the effects of lensing [12], or in the properties of accretion 
phenomena [131 [HI [IS] • 

One of us recently suggested that observations of a hypothetical class of stars 
orbiting very close to the galactic center black hole could provide tests of the so- 
called "no-hair" theorem of general relativity [16]. Specifically, measurements of 
the precessions of the orbital planes of a number of stars with an accuracy of 10 
microarcseconds (/ias) per year could determine both the angular momentum J and the 
quadrupole moment Q of the black hole, and thereby test the constraint Q = —,P/Mc^ 
imposed by the Kerr solution of general relativity. Detection of such stars and achieving 
the required astrometric accuracy are goals of the next-generation of near-infrared, 
adaptive optics interferometry being designed and built by a number of groups [171 [18] . 

However, in assessing the feasibility of such strong-field GR tests, one must 
inevitably address potential complications, notably the perturbing effect of the other 
stars that may also reside in a cluster close to the black hole. Using iV-body simulations, 
Merritt et al ([19], hereafter referred to as MAMW) showed that for a range of possible 
stellar and stellar-mass black hole distributions within the central few milliparsecs (mpc) 
of the black hole, there could exist stars in eccentric orbits with semi-major axes less than 
0.2 milliparsecs for which the orbital-plane precessions induced by the stars and black 
holes would not exceed the relativistic precessions. These conclusions were gleaned from 
thousands of simulations of clusters ranging from seven to 180 stars and stellar mass 
black holes orbiting a 4 x 10^ Mq maximally rotating black hole, taking into account the 
long-term evolution of the system as influenced by close stellar encounters, dynamical 
relaxation effects, and capture of stars by the black hole. 

In this paper we study the extent to which the conclusions of MAMW can be 
understood, at least within an order of magnitude, using analytic orbit perturbation 
theory. We calculate the average change in the orientation of the orbital plane of a 
given "target" star orbiting the massive black hole, as determined by its inclination 
and ascending node angles i and f2, induced by the Newtonian gravitational attraction 
of a distant third star (which could be either inside or outside the target star's orbit). 
The perturbing accelerations are expanded in terms of multipoles through £ = 3. We 
then calculate the root-mean-square variation of each orbit element, averaged over all 
possible orientations of the perturbing star's orbit, and averaged over a distribution of 
orbits in semi-major axis and eccentricity, arguing that this will give an estimate of 
the "noise" induced by the graininess of the otherwise spherically symmetric perturbing 
environment. 
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Figure 1. R.m.s. precession dO/dt ~ {{P) + sin^ for a target star with 

e = 0.95 plotted against semi-major axis, for three models with 7 = 2, /? = 0, 
R — 1. Mi, denotes the total mass within one mpc. Shown (blue in color version) are 
results from Integration I (dashed curves), Integration II (solid curves) and Integration 
III (dot-dash curves). Also shown are the amplitudes of frame-dragging (black in 
color version) and quadrupole (red in color version) relativistic precessions for the 
corresponding star, assuming a maximally rotating black hole. Wide line (orange 
in color version) denotes the precession corresponding to an observed astrometric 
displacement of 10 /^as/yr. 

Figure 1 shows the results for three stellar distribution models, among the set of 
models studied by MAMW. They have number densities that vary as 1/a''', where a' is 
the semi-major axis of the perturbing star, and have eccentricity distributions that vary 
as (1 — e'^)~^; in Fig. 1, 7 = 2 and (3 = 0, corresponding to a distribution with isotropic 
velocity dispersion. They have an equal number of IMq stars and IOMq black holes. 
The label M^, chosen to parallel the notation of [I9], denotes the total mass within 
one mpc of the black hole; the three cases correspond to a total number of perturbing 
bodies within a radius of four mpc of 7, 21 and 72, respectively. The target star has 
eccentricity e = 0.95, and its semi-major axis a ranges from 0.1 to 2 mpc. Plotted is 
the rate of precession of the vector perpendicular to the orbital plane, d9/dt, observed 
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at the source, in arcminutes per year, calculated using three ways of carrying out the 
integrals over the stellar distribution. The dashed line denotes an integration (I) in 
which all perturbing stars are assumed to be sufficiently far from the target star at all 
times that their pericenters are outside its apocenter or that their apocenters are inside 
its pericenter. The solid line denotes an integration (II) in which closer encounters are 
permitted, limited by demanding that all perturbing stars be on orbits such that the 
higher £ contributions to dO/dt be at worst comparable to the contribution at lowest 
order in i. The dot-dashed line denotes an integration (III) which uses a fitting formula 
that interpolates between the extreme limits of a perturbing star well outside the target 
star, and a perturbing star well inside the target star; in this case the integration is over 
the entire stellar distribution. The orange band in each panel denotes the value of dO/dt 
corresponding to an astrometric precession rate dO/dt of 10 /las per year as seen from 
Earth, given by 

{dd / dt) source _ 1.3 ((i9 /(it) Earth /^v 

(arcmin/yr) d (lOyuas/yr) ' 
where a is the semi-major axis in units of mpc; we use 8 kiloparsecs as the distance to 
the galactic center. 

Also plotted are the amplitudes of the frame-dragging and quadrupole precessions 
for a Kerr black hole, given by [T6] 



47r 



GM 



3/2 



c^a{l — 

0.769(1 - e2)~^/^xo"^ arcminyr"^ , (2) 
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GM 



c^a(l — e^) 

^ 7.97 X 10^^(1 - e^y^x^a-^^^Tcmmyi-^ , (3) 

where P = 2iT{a^ /GMY^'^ is the orbital period, and where x = Jc/GM"^ is the dimen- 
sionless Kerr spin parameter, set equal to its maximum value of unity in Fig. 1. 

Because Integration I keeps the stars far from the target star, the precessions are 
small. By contrast, the fitting formula of Integration III is large for very close encounters, 
so not surprisingly, the precessions from that method are large. Integration II gives 
results intermediate between the two. Interestingly, the spread between these methods 
is roughly consistent with the spread between individual precessions obtained in the 
A^-body simulations of MAMW. This can be seen in the top panel of MAMW, Fig. 7, 
which corresponds to the middle panel of Fig. 1 (to properly compare the two figures, 
one must translate between dO/dt and dQ/dt). It can also been in the bottom panel 
of MAMW Fig. 5, where the points labelled by x indicate the mean precessions in the 
absence of black hole spin, for the same three stellar distributions as are shown in Fig. 
1. Thus we regard our three integration methods as giving a reasonable estimate of the 
range of stellar perturbations. 

Comparing the three stellar distributions shown in Fig. 1, we see that the effects 
vary roughly as N^^"^ oc M^^, as expected, from the nature of our r.m.s. calculation. 
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We consider eight different stellar distribution models, and for seven of them, 
consider models with equal numbers of stars and black holes, and models with only 
stars, totalling 15 models. In all but one case, the precessions are generally smaller 
than the ones shown in Fig. 1, and that case is a centrally condensed model with a 
non-isotropic velocity dispersion leading to a preponderance of highly eccentric orbits. 
We conclude that, for a target star in a very eccentric orbit with a < 0.2 mpc, there is a 
reasonable possibility of seeing relativistic frame- dragging and quadrupole effects above 
the level of 10 /ias/yr without undue interference from stellar perturbations. 

The remainder of this paper gives the details underlying these results. In Sec. [21 we 
describe the basic orbital perturbation theory leading to the orbit-averaged variations in 
the orbital elements of a target star. In Sec. [3] we calculate the r.m.s. variations of the 
elements i and Q by averaging over distributions of perturbing stars. Section H] gives the 
numerical results and compares them with those of MAMW. Concluding remarks are 
made in Sec. [51 Appendix A lists the higher-order contributions to the r.m.s. variations. 
Appendix B derives the minimum distance from the black hole reached by a body that 
avoids either tidal disruption or direct capture, and Appendix C shows that the effects 
of tidal deformations on the orbital planes of stellar orbits are negligible. 



2. Orbital perturbation theory at the galactic center 
2.1. Basic equations 

In Newtonian theory, the acceleration a = ai — a2 of a target star with mass rui relative 
to a massive black hole (MBH) with mass m2 in the presence of a perturbing star with 
mass rrts is given by 

^ _ Gm2ri2 GmsVis Gmiri2 ^ Cmgras 

/v»3 i>i3 /v*3 />^3 ^ ^ ' 

'12 '13 '12 '23 

where G is Newton's constant, Vab = — Vb, and rat = \rab\- 

For a perturbing star inside the orbit of the target star ("internal" star), with 
^"32 ^ ^^12 5 Eq- Q can be expanded as 

• G{mi +7712-1- ?7i3)r* Gm^R 



r3 R^ 



i=i ' ^ ^ 



where r = ri2 and R = r23; the capitalized superscripts denote multi- indices, so that 
R^ = R^R^ . . . R^'., and similarly for the partial derivatives; (. . .) denotes a symmetric 
trace- free product. 

For a perturbing star outside the orbit of the target star ("external" star), with 
^"12 ^ ^"23, the expansion takes the form 



G{mi + m2)r' 



e=i ' ^ 



(6) 
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Because rrii <^ m2 and ^ m2, and because we are only concerned in what follows 
with orbital plane effects, we can replace both mi + m2 and mi + m2 + 777,3 with a single 
M, effectively the mass of the MBH. 

Establishing a reference XY plane and a reference Z direction, one defines the 
standard "osculating" orbital elements. The inclination i relative to the reference plane, 
and the angle of ascending node fl between the X axis and the line where the orbital 
and reference planes intersect, fix the orientation of the orbital plane in space. The 
semi-major axis a, eccentricity e, and pericenter angle u determine the orbit in the 
orbital plane. The true anomaly f = (p — u is measured in the orbital plane from the 
pericenter to the location of the body. It is also useful to define an auxiliary angle of 
pericenter w = u + Q cos i which represents a kind of angle measured from the reference 
X-direction, rather than from the nodal line. 

The unit vector n pointing from the MBH to the target star, and the orthogonal 
unit vectors A and h are given by 

n = ex [cos {u + /) cos Q — sin {u + /) sin Q cos i] 

+ ey [cos {u + /) sin Q + sin [u + /) cos Q cos i] 

+ e^[sin {u + /) sini] , 
\ = — ex [sin {u + f) cos Q + cos {u + f) sin Q cos i] 

— ey [sin [u + /) sin Q — cos {u + /) cos Q cos i] 

+ ez[cos {u + /) sini] , 
h = sin z (ex sin Q — ey cos Q) + ez cos i , (7) 

where h is normal to the orbital plane. We also have the osculating orbit definitions 
r = + ecos/), h = \r X v\ = {GMpy^"^, d(f)/dt = h/r'^, and p = a(l — e^) for the 
target star, and R = p'/{l + e'cosF), h' = \KxV\ = {GMp'f^, d(p'/dt = h'/r'^, and 
p' = a'{l — e'^) for the perturbing star, along with its orbital elements i', Q' and u'. 

Defining the perturbing acceleration to be everything in Eqs. ([5]) and ([6]) except 
the leading acceleration —GMr/r^, and defining TZ, S, and W to be the components 
of the perturbing acceleration along n. A, and h respectively, the Lagrange planetary 
equations for variations with time of the target star's orbital elements are given by (see, 
eg. Sec. 7.3 of [20]), 

n7 T 

— =W-cos(u; + /), (10) 
dVL r 

— = W-sin(a; + /)/sinz, (11) 

— =-n—cosf + S — ^sm/. (12) 
dt en en 
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We will work in first-order perturbation theory, whereby we express TZ, S and W in 
terms of osculating orbit variables, set the orbit elements equal to their constant initial 
values in the right-hand side of Eqs. ([H]) - ( fT2l) . and then integrate with respect to time. 

2.2. Time averaged variations in orbit elements 

We want to use the above equations to calculate the time averaged rates of change of 
the orbital elements of the target star, given by dx/dt = T~^ {dx/dt)dt, where T is 
the longest relevant timescale, and x is the element in question. For an external star, T 
would be the orbital period of the perturbing star, while for an internal star, it would 
be the period of the target star. Assuming that the shorter period Ps is much shorter 
than the longer period Pl in each case, then it is straightforward to show that, modulo 
corrections of order of Ps/Pl, 

llx 1 , ,2,3 /"^^ /"^^ 1 ,^ 

where P is the orbital period of the target star, / and F are the true anomalies of the 
target and perturbing stars' orbits, respectively, e' is the eccentricity of the perturbing 
star's orbit, and dx/df = {r'^/h)dx/dt, valid to first order in perturbation theory. 

By way of illustration, we show here the time-averaged changes of orbital elements 
for the i = 1 term induced by an external star [Eq. (|5])], for the special case i' = and 
n' = 0: 



0, (14) 



da 

1^ 15^ e(l-e'2)3/2 . 

-JTt (1-6^)5/2 smcucoscusm t, (15) 

'dt 15 (1 - e'^)^/^ 

— = -Be^t-r-, ^TTT^e^ sinwcoswsinicosi , (16) 

dt 4 (1 — e^j'/^ 

In 3„ (l-e'2)3/2 , 2 r 2 2 N 

^ = - 1^.., _ (1 + 4e^ - 5e^ cos^ u) cos z , (17) 

3 (1-6^^)3/2 2 o^, 2.. 2 2.^ 

— = -^extT:; r-T73-(5cos w - 3 + 5 cos ism a; - cos ?) , (18) 

dt 4 (1 — e^)^/^ 

where i?ext = i'^ir / P){ms/M){p/p')^ . For arbitrary orientations i' and Q' the expressions 
are much more complicated. We have also found the analogous expressions for the i = 2 
and i = 3 terms. These are smaller than the i = 1 results by factors of p/p' and {p/p')"^, 
respectively. 

For an internal star, the i = 1 term of Eq. (jS]) contributes no time-averaged variation 
of any of the elements. The £ = 2 contributions scale as i?int = (2vr/P)(m3/M)(p'/p)^, 
while the i = 3 contributions are smaller by a factor of p' /p. Again, the general 
expressions are long, so we will not display them here. 

Since the orbital energy of the target star is proportional to 1/a, Eq. f|T^ simply 
reflects the absence of a secular energy exchange mechanism between the target and 
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perturbing stars at first order in the perturbations. As a side remark, Eqs. f llSp and 
f|T6|) together imply that (1 — e^)^/^cosz is a constant, so that a decreasing inchnation 
produces an increasing eccentricity; in planetary dynamics this is known as the Kozai 
mechanism [21] . 

3. Perturbations by a distribution of stars 

3.1. Average over orientations of perturbing stellar orbits 

With the time-averaged changes in the orbital elements due to one perturbing star in 
hand, we now turn to the changes caused by a distribution of perturbing stars. We will 
assume a cluster of stars whose orbital orientations (z', Q', u') are randomly distributed. 
We will discuss the distributions in a' and e' later. The "orientation-average" of a 
function will be defined by 

PTT /•2tt /•2tt 

(F) = — - / sin i' di' / d^l' / dcu' F(i' ,uj') . 
Stt^ Jo Jo Jo 

We then find that {dx/dt) = for all four orbit elements e, i, Q and vu, for both internal 
and external stars. The reason is easy to understand: the averaging process is equivalent 
to smearing the perturbing stars' mass over a concentric set of spherically symmetric 
shells. The target star will thus be moving in what amounts to a spherically symmetric, 
1/r potential and its orbit elements will therefore be constant. 

But for a finite number of stars, the potential will not be perfectly spherically 
symmetric, even if the orientations are randomly distributed. It is the effect of this 
discreteness that we wish to estimate. We do this by calculating the root-mean-square 
(r.m.s.) angular average [{{dx/dt)'^)]^^'^. This will give an estimate of the "noise" induced 
in the orbital motion of the target star by the surrounding matter. We will then compare 
this noise with the relativistic effects that we wish to measure. 

We will focus on the quantity 

{(dh/dtf) = {(dijdif) + sin^ i{(dn/dif) , (19) 

rather than the individual elements i and Q; this quantity represents the r.m.s change 
in the direction of h, the normal to the orbital plane. The leading contributions from 
internal and external stars are given by 

{(dhJdtrU = ^BfJ-^^^^^^ , (20) 

{(dhjdt)\.t = ^^cxt|^^^ (l + 3e^ + ye^) . (21) 

Appendix A lists the separate r.m.s. orientation averages for di/dt and dVt/dt for internal 
and external stars, and for all £ < 3. For future use, we define the angular r.m.s. rate 
of change of the orbital orientation by dO/dt = {(dh/dty)^^'^. 
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3.2. Average over size and shape of perturbing stellar orbits 

We now integrate over the semi-major axis a' and eccentricity e' of the perturbing 
stars. We will use a distribution function of the form N'g{a')h{e'^)da'de''^, where is 
normalization factor, set by the condition Af = N/X, where is the total number of 
stars in the distribution, and 

h{e'^)de'^ [ g{a')da', (22) 



where the limits of integration will be determined by the limiting orbital elements for 
those stars. Following MAMW, we will consider a range of parametrized models for 
the dependences g{a') and /i(e'^), and will consider clusters that contain both stars and 
stellar-mass black holes. 

The variables a' and e' will be constrained by a number of considerations. The 
minimum pericenter distance r^^^ for any body will be given by the tidal-disruption 
radius for a star, and the capture radius for a black hole. This will therefore give the 
bound 

a'(l - e') > r^in . (23) 
For rmin we will use the estimates 

,.-«4xlO-'(-) [j^j^J mpc, 

hh 8GM , / M \ 

— -1.5X10- (j^^355^j mpc. (24) 

These are derived in Appendix B. 

However our analytic formulae for the r.m.s. orientation-averaged variations are 
valid only in the limits p' /p ^ 1 or p/p' ^ 1 for internal and external stars, respectively. 
But since our target star is embedded inside the cluster of stars, there may well be 
perturbing stars that do not satisfy either constraint. On the other hand an encounter 
between the target star and another star that is too close could perturb the orbit so 
strongly that it will be unsuitable for any kind of relativity test. Because we are looking 
only for an estimate of the statistical noise induced by the cloud of stars, we will try 
three approaches in order to capture the range of perturbations induced by the cluster. 

Integration I. Because Eqs. fl20|) and fl2T]) are valid only in the extreme limits where 
the perturbing star is always far from the target star (so that the higher-order terms 
are suitably small), we cut out of the stellar distribution any stars that violate this 
constraint. This yields the following conditions on the allowed orbital elements of the 
perturbing stars: (i) for an internal star, we demand that r[^j^^ = a'(l + e') of the 
perturbing star be less than r^in = a(l — e) of the target star; (ii) for an external star, 
we demand that r^j^^ = a'(l — e') of the perturbing star be greater than rmax = fl(l + e) 
of the target star. 

For an internal star, we thus have the two conditions, 

a'{l-e')>r^in, a'{l + e') < a{l - e) . (25) 
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The maximum values of e' and a' allowed under these conditions are 

/ _ ~ ^) ~ ^min I _ 1 ~ 6 , , 

emax,int - _ ^ ^^.^ , «max,int " « ^ + " l^^J 

For an external star, we have the two conditions 

a'(l - e') > a(l + e) , a' < a^nax , (27) 

where amax is the outer boundary of the cluster, chosen to be large enough that the 
effects of stars beyond this boundary are assumed to be negligible. Following MAMW, 
we choose amax = 4 mpc. The maximum e' and minimum a' allowed are thus 
/ ^1 a{l + e) , ^ 1 + e 

^max,cxt 5 "min,ext " -i / ■ 

Thus the average of a function J^{a', e') over this distribution will be given by 

(^) = N{Ji + J2) , (29) 



where 



l2 , 
max,int /" max,int 

Ji(J^) = / h{e'^)de'^ g{a')J^{a',e')da' , 

Jo Jr^in/a-e') 



'rmin/(l-e') 

J2{J^)= I h{e'^)de'^ g{a')T{a',e')da' . (30) 



,/2 

' max.cxt 



However, instead of substituting A/" = N/X, we substitute 

Ar = iV/(Xi+J2), (31) 

where 

pe'^ ■ . pa' 

I max, lilt / '^max.mt 

Ji = / h{e'^)de'^ / g{a')da' , 

Jo Amin/a-e') 

max.Gxt /•Iniax 

X2 = / h{e'^)de'^ / g{a')da' . (32) 

>^«miii.cxt 

This amounts to assuming that all stars in the cluster happen to have orbit elements 
that satisfy our constraint. Thus the average of the function J^{a', e') will be given by 

CF) = Nllin±Mn . (33) 

Note that if J" = 1, we get (J") = N. 

In our simple model, we are treating the stars and black holes as independent 
distributions, so the mean value of J-" can be written as a sum over the two normalized 
distributions, 

{T) = {T)s+{T)b, (34) 

where the only difference between the integrals for the distributions is the value of rmin, 
which affects only the integrals Ji and Xi, and the number of particles, Ns for stars. 
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and Nb for black holes, with = A'"^ + A^^; for later use, we define Nb/Ns = R. Hence 
we obtain 

^ Ns JisjJ^s) + J2{J^s) NbJibI^bHM^ . . 

^ ' N X15+X2 N X1B+X2 ■ ^ ' 

For the r.m.s. variations in dh/dt, we include all the higher-order terms shown in 
Appendix A. 

Integration II. Taking the ratio of the higher i contributions to the orbit element 
variations to the leading i contribution (see Appendix A) reveals that the parameter 
controlling the relative size of the higher-order terms is the ratio a' /a{l — e^) for internal 
stars, and a/a'{l — e'^) for external stars. Requiring each of these ratios in turn to be 
less than one, we repeat the integrals, but with new limits of integration given by 

^max.int ~ ^ ~ ^"111111/^1(1 — 6 ) , '^max.int ~ '^(^ ^ ^ ) ' 

e'max.ext = (1 " a/a^ax)^^^ , a'min.ext =«/(!- e'^) • (36) 

This condition permits closer encounters than the condition imposed in Integration I. 
Here as well, we include all higher-order contributions to the r.m.s. variations. 

Integration III. In an attempt to include even closer encounters between the target 
star and cluster stars, we adopt a fitting formula for the r.m.s. perturbations of the 
orbital plane that interpolates between the two limits of very distant internal and very 
distant external stars. A simple formula that achieves this is given by 

~ {{dh/dtY)7i + {(dh/dty)-^, ' 

where we use only the lowest-order contributions to the r.m.s. variations, given by Eqs. 
f l20|) and f l2T|) . In this case the average over the distributions becomes 

NsM^^NbM^ .... 

where the integrals now take the form 

J(J^) = / h{e''^)de'^ g{a')J^{a\e')da' , (39) 

with X = J(l), thereby including the full distribution of stars. 



4. Numerical results 



In order to compare our analytic estimates with the results of the N-body simulations 
of MAMW, we will adopt as far as possible the same model assumptions. We 
parametrize the distribution functions g{a') and h{e'^) according to g{a') = a'^~'^, and 
h{e''^) = (1 — e'^)~^, where 7 ranges from to 2, and /3 ranges from -1 to 0.5. The values 
(7, /3) = (2, 0) correspond to a mass segregated distribution with isotropic velocity 
dispersion. We will chose a'^^^^ = 4 mpc, arguing that the perturbing effect of the 
cluster outside this radius is negligible by virtue of the increasing distance from the 
target star and the more effective "spherical symmetry" of the mass distribution. We 
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Model 7 /3 7^ N 



Model 7 13 n N 



1 -1 10 159 

2 -1 10 1 29 

3 1 -1 10 119 

4 1 -1 10 1 21 

5 1 30 209 

6 1 30 1 43 

7 2 -1 30 119 

8 2 -1 30 1 21 



9 2 10 1 7 

10 2 30 119 

11 2 30 1 21 

12 2 100 400 

13 2 100 1 72 

14 2 0.5 100 400 

15 2 0.5 100 1 72 



Table 1. Parameters of the distributions 



will assume that the cluster contains stars of mass IMq and black holes of mass lOM©, 
and will consider values of the ratio of the number of black holes to the number of stars 
to be i? = and R = 1 (MAMW also consider the ratio i? = 0.1). The main difference 
between stars and black holes in our integrals is the factor so there will simply be a 
relative factor of 100 between the black hole contribution and the stellar contribution, 
apart from the small effect of the difference in rmin between stars and black holes. 

Of the 22 stellar distribution models listed in Table I of MAMW, we consider only 
the 15 models with either i? = or i? = 1; these are listed in Tabled! While denotes 
the total number of objects within 4 mpc, the parameter denotes the approximate 
total mass within 1 mpc, and gives an idea of the perturbing environment around a 
close-in target star. 

Figure 1 shows the results for the three models 9, 11 and 12, as discussed in Sec. 
[U To illustrate the differences between different models of the stellar distribution. Fig. 
2 shows the predicted precessions for a target star at 0.1 mpc with e = 0.95, for all 15 
model distributions. The crosses and the error bars indicate the range of results from 
the three integration models. Models with 7 = or 1 generally give smaller precessions 
than those with 7 = 2. The latter models are more centrally condensed, and lead to 
larger perturbations of a close-in target star. For the same value of (7, /3, M^), models 
with equal numbers of stars and black holes {R = 1) lead to larger perturbations than 
those with pure stars (i? = 0); the former models are more "grainy" (smaller A^), and 

1 /2 1/2 

so the effects are larger by roughly N^^q/N^I^^. Models 14 and 15 (/3 = 0.5) have an 
excess of stars in highly eccentric orbits, thus leading to larger precessions. 

We have tested the reliability of these estimates by carrying out a number of checks. 
For Integrations I and II, we calculated the precessions first using only the lowest- 
order terms [Eqs. fHOjl . f HTj) . dS]), fH5|) ]. then using those plus the first-order terms 
[Eqs. (USD, g3]), (UnD, gZ])], and finally adding the second-order term [Eqs. (iHD, (l49l)]. 
Table [2] shows, for three values of a and two values of e the percentage change in the 
precession estimates for model distribution 11 caused by adding higher-order terms. 
Not surprisingly Integration I suffers the smallest change, since it forces perturbing 
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Figure 2. R.m.s. precession dO/dt = ((/i^))^/^ for target star with e = 0.95 and 
a = 0.1 mpc for 15 stellar distribution models. Symbol x denotes estimates from 
Integration 11, and error bars indicate the range of estimates from Integrations I and 
III. Amplitudes of frame-dragging, quadrupole, and astrometrie displacement of 10 
/ias/yr are shown as in Fig. 1. 



stars to be always relatively far from the target star, where the higher-order terms are 
relatively smaller. Integration II suffers larger changes, as much as a factor of 2.5 for 
large eccentricities, but still within an order of magnitude. 

We also checked the fitting formula used in Integration III by carrying out a 
calculation of the r.m.s. precession of a target star in an eccentric orbit by a star in a 
circular orbit of the same semi-major axis. The equations of motion can be formulated 
exactly, and the two unperturbed orbits have the same period, allowing a single timescale 
to be extracted. We carried out the time averages of the perturbation equations for 
di/dt and dfl/dt numerically, and then carried out the averages of {dh/dtY over i' and 
VL' numerically. Finally we did a very coarse average over the phase of the circular orbit, 
in order to avoid the singular case where the two stars actually collide! The result was 
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5 
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1 n n 


42.1 


55.8 




0.95 


4.0 


12.2 


49.7 


257.8 


0.5 


0.5 


7.2 


9.4 


41.9 


54.3 




0.95 


3.8 


12.2 


50.3 


255.7 


1.0 


0.5 


6.4 


8.5 


41.3 


53.1 




0.95 


2.8 


12.7 


48.9 


251.7 



Table 2. Percentage change in dO/dt from adding higher order terms in Integrations 
I and II. 

that, for large eccentricity e, the fitting formula overestimates dO/dt by about a factor 
of 2. 

5. Conclusions 

We have used analytic orbital perturbation theory to investigate the rate of precession of 
the orbital plane of a target star orbiting the galactic center black hole SgrA* induced 
by perturbations due to other stars in the central cluster. We found that, although 
the results have a wide spread, they compare well with the distribution of precessions 
obtained using A^-body simulations. One feature not included in our analysis is the fact 
that orbital planes in a real cluster are not randomly distributed, but become somewhat 
correlated over the long-term evolution of the cluster. Whether these correlations are 
large enough to have a significant effect on our estimates is an open question. Within 
our assumptions, however, we find a range of possible models for the cluster of objects 
within the central 4 mpc of the black hole in which it may still be possible to detect 
relativistic precessions of the orbital planes at the 10 yuas/yr level. 
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Appendix A: Compendium of r.m.s. orbital perturbations 

Here we list the r.m.s. perturbations in i and VL for different values of I for both internal 
and external stars. It turns out that cross terms between different I values vanish. 
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Internal: Lowest order {£ = 2) 



,A , _ 3 2 1 + Se'^ + 21e'' 

{[-n) /int - ^Ant — n — — ' y^'^) 



'dt' '™ 80 (l-e'2)4 
-m _ 3 , 1 + 3e-^ + 21e-^ 1 
^^rft^^'"* 80 (l-e'2)4 sin^z' ^ ' 

Internal: First order {£ = 3) 

((^Vt (l-e-^)e + 12 cos-.), (42) 

/S2^ _J[Lr2 /p^y e^e'\6 + 9e^- + 34e^^) (5 + 12 sin^ uj) 
dt ^ ~ 7168 [ pj (1 - e'2)6 sin- z ' ^ 

External: Lowest Order {£ = 1) 

((|^. ^^i^el.|^(C^. + Acos-.), (44) 
dl] 2. _ 3 2 (l-e--)Mgi + Asm-c^) 

External: First order {£ = 2) 

.,^,2, 225 ^2 e-(l -e'2)3^ 

((-r) )cxt = B^^ii — ] -T ^T7^('^2 + ^2Cos-w) , (46) 

^^dt' ' 3584 VpV (l-e2)9 ^ ^ ^ ^' ^ ^ 

5; 2, _^r2 /^pA' e-(l-e-T(C2 + I^2sin-a;) 
rft ^ " 3584^^* U' / (1 - sin- z ' ^^^^ 

External: Second order {£ = 3) 



where 



dt' '"^ 4096 \p' J (l-e2)ii 

X (C3 + D3 cos- w) , (48) 

45 2 /^pA' (1-e--)^ A , o /2 , 7 .-2. 

((— ;- )cxt = -Do^t — 7 1 + 3e + -e sm 2 

dt ' ' 4096 Vp' / (1 - ) V 27 

X (Cs + Z^gsin-w). (49) 

Ci = (1 - e-)- , = 5e-(2 + 3e-) , 

C2 = 5(1 - e-)- , 1^2 = (4 + 3e-)(3 + lie-) , 

C3 = (1 - e-)-(2 + 3e- + 44e^) , = 21e-(2 + e-)(l + 5e- + 8e^) . (50) 

Appendix B: Minimum distance for a stellar or black- hole orbit 

A star that approaches too close to the black hole will be tidally disrupted and be 
removed from the stellar distribution. An estimate of this distance is given by the 
"Roche radius", r^oche ~ R{2M/mY^^ , where R is the radius of the star, and M and 
m are the black-hole and stellar masses, respectively. For a solar-type star, the radius 
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R may be estimated using the empirical formula R ~ RQ{m/mQ)^-^. Thus we obtain 
''"min ~ RQ{m / mo)^-'^'^ (2M / m^Y^^ . Putting in numbers gives the first of Eqs. flM|) . 

A stellar-mass black hole will not be tidally disrupted, but can be captured directly 
if its energy and angular momentum are such that there will be no turning point 
in its radial motion. For equatorial orbits in the Kerr geometry (in Boyer-Lindquist 
coordinates), the equation of radial motion has the form [dr/drY = E"^ — V{r), where 
r is cx proper time, E is the energy per unit mc^ of the body, and 

V{r) = 1 + - + ^ — ' (51) 

where M = GM/c^, a = J/Mc, (3 = Ll — a?E'^, and a = — aE, where J is the 
angular momentum of the central black hole and is the angular momentum per unit 
mc of the orbiting black hole. The critical angular momentum for capture is given by 
that value such that the turning point occurs at the unstable peak of V{r). Since the 
orbiting stars and black holes are in non-relativistic orbits, we can set ~ 1. Under 
these conditions, it is straightforward to show that 



(L,), = ±2M (^1 + V 1 T a/Mj , (52) 

where the upper (lower) sign corresponds to prograde (retrograde) orbits. For a/M = 1, 
the critical angular momenta are 2M and — 2(1 + \/2)M. Converting to the language 
of orbital elements, where = m^GMa{l — e^), we find in the large e limit, 
Ll ^ 2m'^GMrp where is the pericenter distance. The result is that 

/ / \ 2 

„bh 



2M + ^iTa/M j . (53) 

This ranges from 2GM/c^ to ll.QGM/c^ for a/M = 1 and is SGM/c"^ for a = 
(Schwarzschild). We adopt the latter value as a suitable estimate; inserting numbers 
gives the second of Eqs. ([24 



Appendix C: Effects of tidal deformations 

Even if stars survive tidal disruption on passing very close to the MBH at pericenter, 
they will be tidally distorted, and these distortions can affect their orbits. However, we 
argue that, for the stellar orbits of interest, these effects are negligible. For example, 
the rate of pericenter advance due to tidal distortions is given by (Eq. (12.31) of [20] ) 



dcu 307r, M 1 + 3eV2 + eV8 , 



dt P m \ a J (1 - e2)5 

where is the so-called "apsidal constant" of the star, a dimensionless measure of how 
centrally condensed it is. Inserting R = RQ{m/mQ)^'^ , we obtain 
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The variations in i and scale in exactly the same way, but are further suppressed by 
the sine of the angle by which the tidal bulge points out of the orbital plane, resulting 
from the rotation of the star coupled with molecular viscosity, leading to a lag between 
the radial direction and the tidal bulge. This angle is expected to be very small. Thus 
we can conclude that, as far as perturbations of the orbital planes are concerned, tidal 
distortions will not be important. 
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